Multiple Meixner-Pollaczek polynomials 
and the six- vertex model 



Martin Bender* Steven Delvauxj Arno B. J. Kuijlaars^ 
February 22, 2011 



Abstract 

We study multiple orthogonal polynomials of Meixner-Pollaczek type with respect to a 
symmetric system of two orthogonality measures. Our main result is that the limiting dis- 
tribution of the zeros of these polynomials is one component of the solution to a constrained 
vector equilibrium problem. We also provide a Rodrigues formula and closed expressions 
for the recurrence coefficients. The proof of the main result follows from a connection with 
the eigenvalues of block Toeplitz matrices, for which we provide some general results of 
independent interest. 

The motivation for this paper is the study of a model in statistical mechanics, the so- 
called six- vertex model with domain wall boundary conditions, in a particular regime known 
as the free fermion line. We show how the multiple Meixner-Pollaczek polynomials arise in 
an inhomogeneous version of this model. 

Keywords: multiple orthogonal polynomial, Meixner-Pollaczek polynomial, recurrence 
relation, block Toeplitz matrix, potential theory, six-vertex model. 

1 Introduction 

In this paper we study a system of polynomials, orthogonal with respect to two different weight 
functions of Meixner-Pollaczek type. Our work is motivated by the analysis of the six-vertex 
model in statistical mechanics introduced in |24| and studied in many papers since then, see e.g. 
[l[2Il[inil21l2Hll23[3H]- Colomo and Pronko [TUIIII] studied the role of orthogonal polynomials 
in the six-vertex model, and in particular the Meixner-Pollaczek polynomials. An object of major 
physical interest is the partition function of this model. A rather complete analysis of the large n 
asymptotics in the homogeneous case is based on the analysis of a Riemann-Hilbert problem for 
orthogonal polynomials [6, 7, 8, 9J. In the inhomogeneous case, the corresponding reasoning leads 
to questions of asymptotics of multiple orthogonal polynomials, as we will discuss in Section 

Here we will focus on a particular regime of the six- vertex model, known as the "free fermion 
line". In the homogeneous case, this case is trivial, the partition function being identically 1. 
Also in the inhomogeneous case, a closed expression for the partition function can be calculated 
explicitly, so the present results do not give any new insights into the original model, but should 
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rather be considered as providing exact and asymptotic information in its own right on the 
associated polynomials, the multiple Meixner-Pollaczek polynomials. 

An important tool in the study of (usual) orthogonal polynomials P„ on the real line, ap- 
propriately scaled in order to have n-dependent weights, is that their limiting zero distribution 
satisfies an equilibrium problem. This equilibrium problem is an important ingredient for the 
steepest descent analysis of the Riemann-Hilbert problem for orthogonal polynomials, thereby 
allowing to obtain strong and uniform asymptotics of the polynomials, see e.g. |15l I16| . 

For the case of multiple orthogonal polynomials, however, no general result is known about 
the existence of an equilibrium problem. The aim of this paper is to obtain such an equilibrium 
problem for multiple Meixner-Pollaczek polynomials. The equilibrium problem will be posed in 
terms of a couple of measures {vi, V2) and it involves both an external source acting on vi and a 
constraint acting on V2. This structure is very similar to the equilibrium problem for the GUE 
with external source model, [S]. 

Our method for obtaining the equilibrium problem is similar to that applied to other systems 
of multiple orthogonal polynomials |18l 126) . but a distinguishing feature is the characterization 
in terms of the eigenvalue distributions of block Toeplitz matrices, rather than usual (scalar) 
Toeplitz matrices. We will provide some general results on this topic that are of independent 
interest. Along the way we also obtain Rodrigues formulae and closed expressions for the recur- 
rence coefRcients for multiple Meixner-Pollaczek polynomials. 

Inspired by the scalar case [151 116| , one might hope that the equilibrium problem presented 
in this paper can be used in the steepest descent analysis of the Riemann-Hilbert problem for 
multiple Meixner-Pollaczek polynomials, thereby obtaining strong and uniform asymptotics for 
these polynomials. This approach, although interesting, will not be carried out in this paper. We 
also hope that our equilibrium problem might serve as an inspiration to obtain similar results 
for the general inhomogeneous six-vertex model, not necessarily on the free fermion line, which 
in turn could serve as the first step in finding the large n asymptotic analysis for this model. 



2 Statement of results 

2.1 Multiple Meixner-Pollaczek polynomials 

Let wi and 'W2 denote two distinct weight functions of Meixner-Pollaczek type; Wj : R — > M+ 
with 

w,{x)^^e^'^-\T{\ + ix)\\ i = l,2, (2.1) 

where ti, ^2 G (~''''/2, 7r/2), ti ^ t2 and A > are fixed parameters and F denotes Euler's gamma 
function [IJ. Note that in (|2.ip the gamma function is evaluated in a complex argument and that 

|r (A + ix)\^ F (A + ix) F (A - ix) , for any X,x eM.. 

Furthermore, for A fixed, |F(A + i.T)|^ ^ ^-'^l^l as |a;| 00, so the restrictions on tj guarantee 
that Wj is exponentially decaying for x — > ±00. 

Lemma 2.1. (Existence, uniqueness, real and interlacing zeros) For any nonnegative integers ki 
and k2, there is a unique monic polynomial PkiM of degree k = ki+k2 satisfying the orthogonality 
conditions 

/oo 
Pkuk2i^)^''^'Wj{^)^^ = O7 for m = 0, . . . , fcj - 1, j = 1,2. 
-00 
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The zeros of these polynomials are real and interlacing, in the sense that each PkiM has 
k distinct real zeros x\^''^'^ < x^'"^'^ < ■ ■ ■ < x^"^'^ such that x^^"^'^ < x^^^^''^^ < and 

ki .k2 ^ ki .k2 — 1 ^ ki .k'l i i ^ ■ ^ ; i 

X A < X A < x,Vi whenever 1 < 7 < fc — 1. 

Lemma 12.11 will be proved in Section 13.11 

In analogy with the case of Meixner-Pollaczek orthogonal polynomials [23] we refer to the 
PkiM multiple Meixner-Pollaczek polynomials; see also [4l Sec. 4.3.3]. These polynomials are 
related to the six- vertex model; see Section For information on other systems of multiple 
orthogonal polynomials in the literature, see e.g. pi IM]. 

In this paper we will derive a Rodrigues type formula for the polynomials Pki.k2i enabling us 
to compute explicitly the following four term recurrence relations, to be proved in Section [31 

Theorem 2.2. (Recurrence relations) Let ti,t2 € (— 7r/2,7r/2) with ti ^ <2- Then, for non- 
negative integers ki and k2, the multiple Meixner-Polleczek polynomials satisfy the recurrence 
relations 



Pi 



fei+i 



- - 4V*4) P>^^M{^) - b'^;%Pk^M-i{^) - cl';XPk,-iM-ii^), (2.2) 



Pi 



h,k2+ii^) - - Pk^Mi^) - bltk,Pk,-iMi^) - 4-*,\Pfe,-i,fe,-i(a;), (2.3) 

t,^t2 (k + ki+2X) , 
t^^t^ _ (fc + 2A - 1) / fci 



where 



h 

ti,t2 _ fci(fc + 2A - l){k + 2A - 2)(tanii - tania) 
^'^i-'^^ ~ Scos^ti ' 

k — ki + k2, and where we set P-i.^j = 0, Pfcj,_i = for any fci, ^2- 

The above theorem is the basis for the main purpose of the paper, namely to study the 
asymptotic zero distribution of the appropriately rescaled multiple Meixner-Pollaczek polynomi- 
als. Fix ti and ^2- For simplicity we choose a particular sequence of indices {ki{k), k2(k))^i 
along which we analyze the large k asymptotics of Pki(k),k2(k)^ form the single sequence 



Qk(x) 



Pk/2,k/2{x) if fc is even, 

P(k+i)/2,(k-i)/2{x) if fc is odd. 



of polynomials. We will show that the zero-distribution of the rescaled polynomials Qk{kx) has 
a weak limit as k goes to infinity, and in order to study this we let n G N and introduce the 
doubly indexed sequence of monic polynomials {Qk.n}k>o defined by 

Qk,n{x) = \Qk{nx). 
From (j2.2p and (j2.3p we obtain an explicit recurrence relation for Qk.m 

xQk,n{x) = Qk+l,n[x) + ak^nQk,n{x) + bk,nQk-l,7i{x) + Ck,nQk-2,n{x) , (2.4) 
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with initial conditions Q-3,„ = Q-2,n = Q-i.n = 0, where 



1 t, t, 3fc + 4A k 
ak,n = -aj.;, 1./2 = — -, tanii + — tani2 

fc(fc + 2A- 1) / 1 1 



_ J_ tl,t2 _ 



fc(fc + 2A- l)(fc 



cos^ii ^ cos2t2 y 

2A - 2) (tanti - tant2) 

cos^ti 



(2.5) 



for k even and 



_ j_ t2-tl 

'^fe." ~ ^"(fe-l)/2,(fe+l)/2 



fc + 1 3fc + 4A - 1 
— - — tan i 1 H tan t; 



- ^-(fe_l)/2,(fe+l)/2 4^ uc...., , 

_ fc + 2A-l / fc + 1 fc-l\ 

"fe." ~ „2V-l)A(fc+l)/2 - 8^2 I^C0s2ti COs2t2y 

_ J_ t,,t, _ (fc-l)(fc + 2A-l)(fc + 2A-2 

~ ^3'^(fc-l)/2,(fc+l)/2 - — - 



COS^ti C0s2f2 / 

!A - l)(fc + 2A - 2) (tant2 - tanii) 



(2.6) 



for k odd. 

Using the recurrence relation (|2.4p . standard considerations show that the zeros of Qn,n are 
the eigenvalues of the 4-diagonal matrix 

ao.n 1 \ 

C2,n &2,n 02,0 



V 



(2.7) 



Cn — l,n ^n — l,n ^n — l,n / 



So the problem of the asymptotic zero-distribution of reduces to finding the eigenvalue 

asymptotics of (I2.7p . If fc and n both tend to infinity in such a way that k/n s for some 
constant s, then the coefficients a^^n have two subsequential limits and for even and odd k 
respectively. Similarly, bk^n and Ck.n have subsequential limits 6^, and 6^, along subsequences 
consisting of even and odd k respectively. Using the identity l/cos^t — 1 + tan^f, these limits 
become 

flg = (3 tan<i + tant2) s/4 

bs= {2 + tan^ti + tan2t2) s^/8 (2.8) 
Cs = (tan<i — tani2)(l + tan2ii)s'^/16, 



(is = (3tan<2 + tanti)s/4 
is = bs 

Cs = (tanf2 — tanii)(l + tan2t2)s^/16. 
If we consider (|2.7I) in 2 x 2 blocks, it is tri-diagonal with blocks 



(2.9) 



A 





1 



A 



(0) 



a2fe,; 



1 



b2k+l,n 0'2k+l,' 



and 



A 



(1) 



C2fe,; 



b7h 



\ C2k 



4 



In view of the limits (I2.5P and (12. 6p . the blocks are slowly varying along the diagonals if n is 
large. In other words, for large n, the matrix (|2.7p has a locally block Toeplitz structure with 
square blocks of size r — 2. This allows its limiting eigenvalue distribution to be obtained from 
a general machinery to which we turn now. 



2.2 Polynomials generated by a general recurrence relation 



/ Qo,nix)\ 




/ Qo,nix)\ 


Ql,n{x) 




Ql,n{x) 


V : / 







In this subsection we will work in the following general setting. Let n be a fixed parameter 
and let (Q k,n{x))'^Q be a sequence of monic polynomials, where Qk,n has degree k and depends 
parametrically on n. Assume that the Qk,n are generated by the recurrence relation 



(2.10) 



where J„ is a semi- infinite matrix with unit lower Hessenberg structure, i.e., the strictly upper 
triangular part of J„ is equal to zero, except for the first superdiagonal, on which all entries 
are 1. We also assume that the lower triangular part of J„ has a finite bandwidth, which is 
independent of n. 

The entries of J„ are assumed to have asymptotically periodic behavior with period r (r > 1). 
More precisely, suppose J„ is partitioned into blocks of size r x r, with one superdiagonal and a 
fixed finite number /3 of subdiagonal non-zero blocks. 



Jn 



(-1) 
0,n 

,(0) 







^l,n ^1 



(-1) 



Am 

0,n 



/3,n 



A 



(/5-2) 
(/3-1) 



(2.11) 



where 



A 



(-1) 

k,7i 





\1 





0/ 



/* 1 



A 



(0) 
k.n 



(2.12) 



where the *'s denote arbitrary constants. 

We will assume that the block entries A^^^^ in (|2.1ip are slowly varying with n, in the sense 
that the limits 

i'^') ~ A^^ (2.13) 



lim 



A 



k,n 



exist for any j = —1,0,...,/? and s > 0. Here the notation lim, 



k/r, 



means that we let both 



k and n tend to infinity, in such a way that the ratio rk/n tends to a limit s > 0. The relation 
(|2.13p is on the level of r x ?■ matrices, with the limit taken entry-wise. If (I2.13P holds then the 
matrix in 1)2. lip is said to have locally block Toeplitz structure, in the spirit of [33] . 
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For fixed s, we collect the limiting matrices in (|2.13p into the following matrix- valued Laurent 
polynomial As{z): 

A,{z) Ai-^h-^ + + • • • + Ai^h^. (2.14) 
We will sometimes refer to As{z) as the symbol. From (|2.12p - (|2.14p it follows that 

/ * 1 0\ 



Asiz) 



* 1 

y(z^"'" + *) * • • • */ 



-Oiz), (2.15) 



where 0{z) denotes all the terms in (12.141) that tend to zero as z — 0. 
Define for each s > 0, the algebraic equation 

fs{z, x) det(A,(z) - xir) = 0, (2.16) 

where denotes the identity matrix of size r. Note that fs depends on two complex variables 
z and X. By expanding the determinant (j2.16|) and using p.lSp . we can write fs as a (scalar) 
Laurent polynomial in z: 

Mz, x) = /_i,,(x)z-i + /o,,(x) + ■ ■ • + /p,,(x)zP, (2.17) 

where /j^s(x), j = —1, 0, . . . ,p, are polynomials in x, with /_i^s(a;) = (— 1)*"^^. We define p in 
(|2.17p as the largest positive integer for which /p.s ^ 0. 
Let us solve (j2.16p for z; this yields p + I roots 

Zj = Zj(x, s), j = l,...,p+l. 

We assume that for each x G C these roots are ordered such that 

< ki(x,s)| < |z2(a;,s)| < ••• < \zp+i{x,s)\. (2.18) 

If cc € C is such that two or more consecutive roots in (I2.18P have the same absolute value, then 
we arbitrarily label them so that (|2.18p is satisfied. 
Define the set 

ri(s) = {x€C \ \ziix,s)\ = \z2ix,s)\}. (2.19) 
In the cases we are interested in, we will have that 

ri(s)cR. (2.20) 

Supposing this to hold, we define a measure /if on ri(s) C M with density 

d^?(x) = i-L f 'A±^ - dx. (2.21) 

' r2'Ki \zi+{x,s) zi^{x,s)) ^ ' 

Here the prime denotes the derivative with respect to x, and zi±{x, s) are the boundary values of 
zi{x,s) obtained from the +- side (upper side) and — side (lower side) respectively of ri(s) C K. 
These boundary values exist for all but a finite number of points. 

As discussed in [17|, the measure /if can be interpreted as the weak limit as n — > cx) of the 
normalized eigenvalue counting measures for the block Toeplitz matrices Tn{As) associated to 
the symbol (pHll . 
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Lemma 2.3. With the above notation, we have 

(a) zi{x, s) — x^"^ + O (x^''^"'") as X ^ oo. 

(b) /if in (|2.2ip is a Borel probability measure on Ti{s). 

Proof. See P7]. □ 

With this notation in place, let us return to the sequence of polynomials {Qk,n)kLo '^^ (|2.10p . 
This sequence is said to have real and interlacing zeros if each Qfc,n has k distinct real zeros 

k,n ^ k,n ^ ^ k.n i i-i , k.n ^ k — l,n ^ k.n i i ^ ■ ^ 7 1 

Xi < X2 < ■ ■ ■ < Xf. such that < < "^j+i whenever 1 < j < fe — 1. 

The next result states that, under certain conditions, the normalized zero-counting measures 
of the polynomials Qn,n have a (weak) limit for n — > oo, which is precisely the average of the 
measures (j2.2ip . Here the average is with respect to the parameter s. 

Theorem 2.4. (Limiting zero distribution of Qn ^i) Let the sequence of polynomials {Qk.n)kLo 
be such that (j2.10p - (|2.13p hold. Assume that {Qk.n)kLo ^'^^ '"^'^^ '^'^'^ interlacing zeros for each 
n, as described above. Also assume that (|2.20p holds for every s > 0. Then as n ^ oo, the 
normalized zero-counting measure p„ of the polynomial Qn,n, 

^-=^E%"' (2.22) 
where Sz denotes the Dirac point mass at z, has the limit 

lim pn — vi I III ds (2.23) 

in the sense of weak convergence of measures, where p,\ is defined in (j2.2ip . 

Theorem 12 .41 will be proven in Section This theorem generalizes a result for the scalar case 
r = 1 by Kuijlaars-Roman Theorem 1.2]; see also [T31[57j. 



2.3 Zeros of multiple Meixner-Pollaczek polynomials 

We apply Theorem 12.41 to the polynomials Qk.n iu Section \m In what follows, we will assume 
the condition of symmetric weights, 

t :— ti ~ —t2\ 

for this case we can characterize the limiting zero distribution in terms of a constrained vec- 
tor equilibrium problem, Theorem 12.61 which is the main result of the paper. The symmetry 
condition in Proposition 12 . 51 is needed to prove that Ti(s) C K; in the general case this may fail. 
The recurrence coefhcients in (|2.8p and (|2.9p then become 

Og = — flg = (tant)s/2, 

bs = Is = (1 + lSinh)s^/A, (2.24) 
Cs — — Cs = tanf(l + tan^<)s'^/8 = agbs. 

We partition the matrix J„ from (12. 7p into blocks as in (|2.1ip with r = 2 and /3 = 1. Then 
the limiting values in (j2.13p exist and are given by 

As -(^1 oj' 
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and 



as 1 

—a. 



-c. 



^(1) 

with Us, bs and Cs as in (|2.24p . The symbol As{z) in (|2.14p now becomes 
so (|2.16p reduces to 

- + + . 0. (2.25) 

This equation has three roots, Zj{x, s), j ~ 1,2,3 which we label in order of increasing modulus: 

< \zi{x, s)\ < \z2{x, s)\ < \z3{x, s)\. (2.26) 

Proposition 2.5. (Limiting zero distribution) Suppose ti = — 12 = i G (0,7r/2). As n — >■ oo, 
the normalized zero-counting measure pn of the multiple Meixner-Pollaczek polynomial Qn.n in 
Section \2.1[ see (|2.22p , converges weakly to the average vi = fif ds of the measures in 

(Ml. 

Proof. This is a consequence of Theorem l2.4l The assumptions in the latter theorem are indeed 
satisfied: the interlacing condition follows from Lemma \TJ\ and the fact that ri(s) C M follows 
from Proposition 15. II □ 

In addition to the set Ti{s) in (j2.19p . we define 

r2(s) = {xeC\ \z2ix, s)\ = \z3{x,s)\}. (2.27) 

We will now characterize the limiting zero distribution vi in terms of a vector equilibrium prob- 
lem from logarithmic potential theory |30l I31| . Recall that for a pair of Borel measures (/i, j^) 
supported in the complex plane, the mixed logarithmic energy of fj, and v is defined as ^31j 



/(/i, log 1^ _ ^1 dM(x) Av{y). 

Theorem 2.6. (Equilibrium problem) Suppose ti — — <2 ~ t £ (0,7r/2). Then the asymptotic 
zero distribution, vi, ofQn^n is the first component of the unique minimizer (t^i,;/2) of the energy 
functional 

Eifi, v) := fi) + I{v, v) - {tt - 2t)\x\ dp{x). 



among all vectors (/x, i^) of positive measures such that supp /i C M, J d/i = 1, suppi/ C 
J dv = 1/2 and v is absolutely continuous with density satisfying 

di'{ix) ^ 2t 

\dx\ ~ TT 

The measures have the properties 

supp 1^1 = [-Ci,Ci] 
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Figure 1: Density of the measure vi for (a) h — tant = 100 and (b) h = tani ~ 0.1. 



and 



supp((T - V2) = iM \ {~ic2, 102), 



where a is the positive, absolutely continuous measure on iR with constant density ^ , and where 
Ci and C2 are positive constants given by 



ci = 



27b^ + 18P ~ 1 + V62TT(96^ + 1)3/2 
3262 



1/2 



(2.28) 



C2 = 



V^TT{9b'^ + l)^/^ - 276^ - 1862 + 1 
3262 



1/2 



(2.29) 



with 



b = tant. 

Furthermore, vi is absolutely continuous with respect to Lebesgue measure and has density 

1 + iw{x) 



dz^i 1 , 
— — = — log 

da; 27r 



X £ [-ci,ci], 



whe 



w{x) = 



1 — iw{x) 

(4 + z(62-l))(4 + z(62 + l)) 
IGzlxl 



and z = z{x) is the complex solution to the algebraic equation 

(4 + z(62 + l))2(4 + z62) ^ , 
64z ~ ^ 



(2.30) 



(2.31) 



(2.32) 



such that lm{w{x)) < 0. 
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Remark 2.7. In general, there is no closed expression for the density of i^i, but in the limiting 
case t — > 0, corresponding to the ordinary Meixner-Pollaczek orthogonal polynomials, (|2.32p can 
be explicitly solved, z — — 4 + 8i|a;|Vl — 2;^, giving \z\ = 4. Thus (I2.3ip becomes 

, , 16 — 16z— IzPz Im(z) r 

16z|a;| 16|zp|a;| 8|a;| 

and the density (|2.32l) takes the form 



d^^2^^"g(^ l-Vl-.^ j^^'-'^^^- ^'-''^ 

For non-zero values of t e (0,7r/2), the density of vi turns out to have the same qualitative 
features as p.BSp : see Figure [T] for some illustrations. 

Theorem 12.61 will be proven in Section [5] The proof uses a general result for block Toeplitz 
matrices [T7], see also [TS], together with some explicit calculations using (|2.25p . 



2.4 Motivation: The six-vertex model with domain wall boundary con- 
ditions 

Multiple Meixner-Pollaczek polynomials appear in the study of the six- vertex model in statistical 
mechanics, as we explain now. 

Consider an N x N square lattice in the plane. A configuration of the six-vertex model is 
an assignment of an orientation to the edges of the lattice in such a way that each vertex is 
surrounded by precisely two incoming and two outgoing edges. See Fig. [2] for a configuration 
with N — 5. The name six-vertex model refers to the fact that the local behavior near each 
vertex is given by six possible edge configurations (see Fig. [3]). 

We consider the six- vertex model with domain wall boundary conditions (DWBC). This 
means that the edges at the top and bottom of the lattice must be directed outwards and those 
at the left and right of the lattice must be directed inwards. 



■X5 



yi y2 ya yA 

Figure 2: A configuration of the six-vertex model with DWBC for N = 5 



~* — Xi 



■X2 



■ Xi 
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1 2 3 4 5 6 

Figure 3: Six types of vertices. 

To each of the N rows of the lattice we associate a parameter Xi € R and similarly to each 
column a parameter ijj € M, i, = 1, . . . , A^. We also fix a positive parameter 7 and we assume 
that \xi — Ujl < 7 for all i and j. We define the weight of the vertex in row i and column j 
according to its type as sin(27) (type 1 or 2), sin(7— (xi —yj)) (type 3 or 4), or sin(7+ (xi — 
(type 5 or 6). Note that the weights are parameterized according to the so-called disordered 
phase convention. 

The weight of a configuration is defined as the product of the weights of all the N x N vertices 
in the configuration. The partition function Zn = Zn{xi, . . . , xat, j/i, . . . , j/at, 7) is defined as the 
sum of the weights of all the consistent configurations of the NxN six-vertex model with DWBC. 
An explicit expression for the partition function in terms of an A'^ x A'^ determinant was found 
by Izergin and Korepin [501 EI], see also p i^ESj . 

The homogeneous case 

In the case where Xi — yj = t for all i, j = 1, . . . ,iV, for some fixed parameter t e (—7,7), the 
Izergin-Korepin formula reduces to 

= Z^(7, t) = Mn(7 + 0Bin(7^^r^ ^^^^ ^2.34) 



n 



IV. 

, n=0 



where the matrix M = {'mi,j)fj^-^ has entries 



with 



x*+J-2e*^u.(x) dx (2.35) 

00 

( ^ sinh(f(7r-27)) 

'^y^) = ■\^x N , 2.36 

smh(i7r) 

see e.g. [6l 1211 [38] . The matrix M is then precisely the moment matrix corresponding to the 
weight function e*^w(x) on the real line. Standard considerations (e.g. [35]) show that detM 
can be expressed in terms of the monic orthogonal polynomials P„(x) defined by 

P„(x) =x"-f ©(x"-!) 

for all n and 

P„(x)x"'e*"«;(x) dx = /i„5„,„ (2.37) 



for all n, m with m < n. In fact, det M is expressed in terms of the numbers h„ in (j2.37p through 
the formula 



detM = Jl hn. (2.38) 



n=0 
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Special choices of parameters lead to known families of orthogonal polynomials. Indeed, 
Colomo and Pronko [lOlfTT] showed that the Continuous Halm, Meixner-Pollaczek and continuous 
Dual Hahn polynomials appear in this way. In more general cases, the expressions (j2.34l) - (j2.38p 
were used to compute the asymptotics of the partition function Zjq for large N in great detail 
by means of the Riemann-Hilbert method [B], see also [71[B]. 



The inhomogeneous case 

The situation in this paper corresponds to the case where 

^^-Vj^i ^-'rj'^V ^ (2-39) 

for some ti ^ t2 and all j = 1, . . . , A^. Following the reasoning in [5^, see also the appendix in 
^12j , one sees that the Izergin-Korepin formula reduces to 

2 ^ Nn(7 + h) sin(7 - ^i)]"^^[sin(7 + ^2) sin(7 - t^T^^ 

^ /"l-l \ /n2-l \ /N-l \ \ ■ ) 



n-! n-! n 



, 11=0 / \ n=0 / \ n=0 / 



where the matrix M — (m^^ 



/ x'+^-2e*i^«;(a;)dx, 1 < i < ni, 
m^j = { ^-.§§ (2.41) 
/ x'^^-'^^-'^e''^'^w{x)Ax, ni + l<i<iV, 
/ — 00 

with w still given by (12.3611 . Thus M is a moment matrix with respect to the system of weight 
functions e*^^w[x) and e*^^w(x) on the real line. 

The inhomogeneous model (|2.40p was studied in |12[ [T5] in connection with the calculation of 
the arctic curve. Leading order asymptotics of the partition function for the case ni = 1 in (|2.40p 
was computed in [T2] for the disordered regime, and in [T3] for the anti-ferroelectric regime. The 
analysis is valid in fact for any ni as long as n\ = o{N) as — 00. 

It turns out that det M can be expressed in terms of monic multiple orthogonal polynomials 
PkiM^^) with respect to the system of weight functions e*i^u;(a;) and e*^^u;(a;). The polynomial 
PkiM^^) defined for any nonnegative integers fci, ^2 by 

and 

poo 

Pk,M{^>""^''"w{^)^x = h^lkMk:^ m = 0,...,fc,, j = l,2. (2.42) 
We now have the following generalization of the formula (|2.38p . 

Proposition 2.8. (Partition function) Let^ > 0, <i,t2 G (^7; 7); ti 7^ ^2 o,nd recall the notation 
w(x) in (|2.36p and h!"^^ in (|2.42p . Then the moment matrix M in (|2.4ip has determinant 

JV-l 

detM=n4f(:L(«)' (2.43) 

n=0 
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1 ifj{n)^i, 
0, otherwise, 



where {j{n))^^Q is any sequence of I's and 2's such that 1 appears ni times and 2 appears 
712 times, and where the associated sequences {ki{n)))^^Q , i = 1,2, are defined recursively by 
ki{0) = and 

h{n) - k,{n - 1) = 
for any 7i = 1 , . . . , iV — 1 and i = 1,2. 

Proof. Wc will give the proof for the particular sequence {j{n))^^^ = (1,...,1,2,...,2); it will 

ni times n2 times 

be straightforward to extend the proof to the more general sequences {j{n))^^Q in the statement 
of the proposition. From the definition (|2.4ip it follows that 

where we define the inner product 

/oo 
f{x)g{x) dx 
-oo 

and where we use the functions 

fr{x) ^ x'e^^^'wix), i = 0,...,ni - 1, 
fm+t{x) = a;'e*='^w(x), i = 0, . . . , ria - 1, 

and 

5„(a;)=a;", n = 0, . . . , TV - 1. 

Define functions i.pi{x) and ipj{^)j *i J = Oj ■ • ■ i ~ li by bi-orthogonalizing the functions fi{x) 
and gj{x) in the following way: 



(p. 



(x) = fi{x) + ^ Ck,ifk{x), 



fc=l 



J-1 



= 9j{x) + hk.jg^{x), 

k=l 

for appropriate coefficients bk,j and Ck,i, subject to the orthogonality relations 



{ip„^j) ^ hiS.^j. (2.44) 

at we can identif 

Then we have that 



It is not hard to see that we can identify = Pki{n).k2in)ix) and hn = h^ki(l) M{n 



N-l 

detM = ((/^,.gj))fj^o i^Vi,'^3))li2o = n 

ri=0 

where the last step follows in a trivial way from (|2.44p . In view of the identifications mentioned 
in the previous paragraph, we then obtain (|2.43p . □ 

It is straightforward to generalize the above reasoning to the case of multiple values of the 
differences in p.39p . In general, one could even allow both Xi and yj to take multiple values and 
then one should deal with multiple orthogonal polynomials of mixed type. We leave the details 
to the interested reader. 
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The free fermion line: Meixner-Pollaczek weights 

The value 7 = 7r/4 corresponds to the so-called free fermion line. As first observed in [ITj, in this 
case the above weight functions are related to the Meixner-Pollaczek weight. Indeed, we then 
have 

e w[x)-e g.^j^(^)-e 2cosh(^)- ^^'^^^ 
We may compare this with the classical Meixner-Pollaczek weight function [23| . 

1 

2^ 



^e^''-'T{x + ix)r{X-ix). (2.46) 



By invoking the identity 

r(z)r(i-z) = -^, 

we see that (j2.46[) for A = 1/2 and t — 2tj reduces to 

J_e«,^p/W2 ^ ix)T(l/2 - ix) = e^*^^ — ^. (2.47) 

27r 2 cosh(7rx) 

Thus the weight functions (j2.45p and (j2.47p are the same up to a scaling of the variable x by a 
factor 4. 

It is an easy job to evaluate the Izergin-Korepin formula for the partition function explicitly 
on the free fermion line, so our results will not lead to new insights in that perspective. They 
should rather be considered as giving exact and asymptotic information on the multiple Meixner- 
Pollaczek polynomials in its own right. 



2.5 Outline of the paper 

The rest of this paper is organized as follows: In Section|3]we consider multiple Meixner-Pollaczek 
polynomials with respect to two general weights and prove Theorem l2.2l In Section|3]we establish 
Theorem 12.41 and Lemma 12.11 which together lead to Proposition 12.51 In Section [S] we prove 
Theorem 12.61 using the theory of eigenvalue asymptotics for banded block Toeplitz matrices. 



3 Proofs of Lemma 12.11 and Theorem 12.2 
3.1 Proof of Lemma 12.11 

In this section we prove Lemma [^?T] on the fact that the multiple Meixner-Pollaczek polynomials 
Pki,k2 exist and are unique, and have real and interlacing zeros. It is well known that the 
corresponding result about orthogonal polynomials with respect to one weight function holds; 
we will rely on a generalization of this fact due to Kershaw, [22]. In our context it amounts to 
the statement that a sufficient condition is that for any non-negative integer k — ki + k2 and 
any polynomials A and B (not both identically zero) of degrees at most ki and ^2 — 1 (or fci — 1 
and fc2) respectively, the function f{x) = A(x)wi{x) + B{x)w2{x) has at most k zeros. Since 

A{x)wi{x) + B{x)w2{x) = -i-e^*^^ |r (A + ix)\^ f A(a;)e2(*i-*^)=^ 4- B{x)) , 

2iT V / 

the conclusion will certainly follow if we can show that, for any real t, g{x) :— 74(x)e^*^ -I- B{x) 
has at most k zeros whenever A and B are polynomials such that degA + degB < fc — 1. (By 
convention, the zero polynomial has degree —1.) This can easily be shown by induction, see e.g. 
[301 P- 138]. □ 
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3.2 Some generalities 

For j = 1, 2, let Wj be integrable real functions on the real line such that the measures 'Wj{x)dx 
have moments of all orders. Suppose that for any non-negative integers fci and fc2 there exists a 
unique monic multiple orthogonal polynomial Pki.k2 with respect to the weights wi, W2, that is, 
a polynomial of degree k = ki + k2 satisfying the orthogonality conditions 

Pk-i^,k2{x)x"^Wj{x) dx = 0, for TO = 0, . . . , fcj — 1, j — 1,2. 

Let jkiM denote the sub-leading coefhcient of PkiM^ that PkiMi^) — x'^ ~\- IkiM^^'"^ + 
0(a;'=-2). For j = 1,2, put 

/OO 
x^^Pk,M{^)w,{x)dx, (3.1) 
-oo 

the first non- vanishing moments. 

We begin by stating a general four term recurrence formula for multiple orthogonal polyno- 
mials on the real line in terms of their sub-leading coefficients and first non- vanishing moments. 
This standard fact is shown in |29) and can also be derived from the Riemann-Hilbert problem 
for multiple orthogonal polynomials, see |35| . 

Proposition 3.1. For any positive integers fci and k2, the multiple orthogonal polynomials satisfy 
the following four term recurrence relation: 



^1,^'2 I k\.k2 

Pfe,„i,fc,„i(z). (3.2) 



""fci-l^fes-l 



By symmetry between the two indices, a corresponding recurrence relation for Pk-^^k2+i{z) is 
obtained by interchanging ki and k2 and superindices. 

3.3 Rodrigues formula for multiple Meixner-Pollaczek polynomials 

From standard results on the ordinary Meixner-Pollaczek polynomials, (see e.g. equations (1.7.2) 
(1.7.4) in [23] with </> = <i + tt/2), we have the orthogonality relation 

Pm,o{x)Pn^o{x)wi{x) dx = "■^(" + ^^)^ ^ (33^ 

COS t]^ 1 

and the recurrence relation 

fifji -|- 2A 1) 

xPnA^) = Pn+i,oix) + {n + A)(tanti)P„,o(a;) + 5- '-Pn-i,o{x). (3.4) 

4cos^ii 

For any real parameter t, define the finite difference operator Lt, acting on functions / : C — > 
C, by the equation 

(Ltf) (x) = e^'fix + i/2) - e-''f{x - i/2). (3.5) 
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Lemma 3.2. For any real ti, t2, the operators Lt^ and L^^ commute. 

Proof. Straightforward calculation. □ 

Lemma 3.3. Let t be a real number and f , g analytic functions in a domain containing the strip 
VL = {z :\ Iin(z)| < 1/2}, and assume that there are positive numbers C and e such that 

\f{z)g{z)e''*'\<Ce-'\^<'^\ (3.6) 

for all z Gfl. Then the following integration by parts formula holds: 

/oo /.oo 
fix) (Ltg) (x)e2*- dx = - (Lof) {x)g{x)e''- dx. (3.7) 
-oo J — oo 

Proof. Using the definition of Li, we can split the integral of the left hand side into two terms 
and shift the contours of integration from the real line to M± «/2 for the first and second terms, 
respectively, by Cauchy's theorem. This gives 

r /(x)e2*- (Ltg) (x) dx = H f{x) (e'^^-+'^g{x + i/2) - e'^^--'^g{x - i/2)) dx 

J — oo J — oo 

{f{u - il2)e^'^g(u) - f{u + i j i)e^''^ g{u)) du 

) 

/OO 
(Lo/) («)e2*"5(«)d«. 
-OO 

□ 

We can now derive a Rodrigues type formula for the multiple Meixner-PoUaczek polynomials, 
which will be the tool to calculate explicit recurrence coefficients. 

Proposition 3.4. Let k\ and k2 be non-negative integers and put k = fci + ^2. Let := 

Lo ■ ■ ■ o L denote the mth iterate of an operator L. Then, for any ti,t2 € (—• 7r/2,7r/2) with 

m times 

ti ^ t2, and any A > 0, the multiple Meixner-Pollaczek polynomial PkiM satisfies the Rodrigues 
formula 

(l^^IL^,1 (|r(A + fc/2 + i ■ )|')) {x)=Ck,MPkukAx)\T{\ + ix)\\ (3.8) 

where 

Ck,M = (-2^)^ {costif'{cost2p. (3.9) 



Proof. Define the function 



fm{x) = |r(A + m/2 + ia;)|', 



for any non- negative integer m. First of all, we note that by the properties r(2;+ 1) = zV{z) and 
r(z) = r(z) of the gamma function, it follows immediately from the definitions that if Rm{x) is 
a polynomial of degree m with leading coefficient a„i, then 

{Lt {Rmfi)) {x) =e'*Rm{x + i/2)T (A + ix) T (A + 1 - ix) 

- e""i?™(x - i/2)r (A + 1 + ia;) r (A - ix) 
= (e" (A - ix) Rra{x + i/2) - e-** (A + ix) Rm{x - i/2)) fo{x) 
=Rm+i{x)fo{x) 
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where Rm+i is a polynomial of degree m + 1 with leading coefficient 

o-m+i — ~2iamCost. (3.10) 
Equivalently, by simply replacing the parameter A by A + n/2, 

{Lt {Rmfn+l)) {x) = i?„,+ i(x)/„(x). (3.11) 

By induction over k, it follows from (13.111) and (j3.10p that 

where Rk is a polynomial of degree k in x with leading coefficient (— 2icost)'^. Therefore, 



{L^tlL^tlfk) [x) = Ck,MPkuk2{x)h{x), 



for some monic fcth degree polynomial PkiM- 

Using this representation we can check the orthogonality conditions. Let ni < fci be a 
non-negative integer. Note that the choices f{z) = |r(c + iz)|^ for any real c > 1/2 and g a 
polynomial satisfy condition (|3.6p of Lemma 13.31 this can be seen from the asymptotics of the 
Gamma function valid as \z\ oo with | Arg(z)| <7r — e, 



r(z) = \/2^z(^)'(l + o(l)) 

(Stirling's formula). By definition of the weight function and applying Lemma 13.31 fci times, we 
get 

Pu,M{^)x"^Mx)dx r x^ (l\IL\1!^ (a;)e2*i-dx 



since Lq acting on non-zero polynomials decreases their degree by one. By Lemma [3?2l the same 
argument applies in checking the orthogonality relations with respect to ■u;2- □ 

3.4 Proof of Theorem [272] 

We will need explicit expressions for the first non- vanishing moments, defined by p.ip . These 
are readily calculated using Proposition 13.41 

Proposition 3.5. The first non-vanishing moments hj:^ of the multiple Meixner-Pollaczek 
polynomial Pki,k2 "'^e given by 

,(1) r(2A + fc)fci!(sin(ti-t2))fe^ 

kiM 22^+'^+fei(costi)'=+'=i+2^(cosi2)'== ■ 

and 

,(2) r(2A-ffc)fc2!(sin(t2-ti))'-^ 

kiM 22^+'=+fe2(cost2)'^+'=^+2^(cosii)'=i ■ ^ ^ ^ 
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Proof. Consider the case j = 1; by Lemma l3T2] the j — 2 case is completely analogous. Reasoning 
as in the proof of the orthogonality relations, and noting that Lq{x^) is a polynomial of degree 
/c — 1 with leading coefficient ik, and that _Lo(e^*) = 2ie^^*sint, we find that 

' / (2zsin(ti -t2))'=^e^-('-*=)e2*=^/,(x)dx 

fci!r(2A + fc)(sin(ii -<2))''" 
- 22A+fe+fci (cos ^1 +2^ (cos ^2 )''^ ' 

Here we made use of the orthogonality relation p. 31) (with m = n = 0) to compute the integral. 

□ 

Proposition 3.6. The sub-leading coefficient jkiM of Pki.k^i^) given by 

(2A + fc-l),, , , , , 

7fci,fc2 = 2 (fcitanti + fc2tant2)- (3.14) 

Proof. Let fc = fci + ^2 be fixed. The polynomial P^j+i^^j — Pfe^^fe^+i clearly of degree k and 
satisfies fcj orthogonality conditions with respect to , for j = 1,2. It is thus a multiple of P^j^fe^, 
and reading off the leading coefficient gives Pk^+iM - PkiM+i = ilki+iM - 7fci,fc2+i)-P'ci.fe2- 
Multiplying this relation by x'^'^ and integrating with respect to wi gives the equation 

~ ^i?.fc2 + l ^ (7fel + l,fe2 - 7fel,fe2 + l)^fei^,fe2' 

which by Proposition 13.51 can be written 

7fei+i./c2 - 7/C1./C2+1 = -^^^^^i^ (tan t2 - tanii). (3.15) 

Identifying coefficients in the recurrence relation p.4p for the ordinary monic orthogonal Meixner- 
PoUaczek polynomials Pki,o with respect to wi, gives 7fei+i,o — 7fei,o + Ofei and so 



ki 



7fci+i,o = -^0.] = -■^-^-^^(2A + fci)tanti. 



Then repeated application of p.l5p leads to the claim, for any ki + k2 — k. □ 
Finally we are ready for the proof of Theorem [ 



Proof of Theorem \2.Sl With the explicit expressions for jki.k2 and hj:^ k2 gi'^^n by Propositions 
and 13.61 this follows from the general recurrence relation (13. 2p . □ 



4 Proof of Theorem 12.41 

In this section we will prove the general Theorem 12.41 on the asymptotic zero distribution of a 
sequence of polynomials Qfe,n generated by a recurrence relation (j2.10p - (j2.13p . The main idea of 
the proof follows the scalar case r = 1 by Kuijlaars-Roman [26, Theorem 1.2], see also [T3', '27]. 
But we will need some nontrivial modifications due to the fact that r may be greater than 1. 

The main tool in the proof is the following result on ratio asymptotics for the Qk,m compare 
with [26, Lemma 2.2]. 
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Lemma 4.1. (Ratio asymptotics) Under the assumptions of Theorem \2.4\ we have that for each 
s > there exists R > so that all zeros of Qk,n belong to [—R,R] whenever k < {s + l)n. 
Moreover, 

-^^-M^=zi(x,.), (4.1) 

k/n^s l^k+r,n(X) 

uniformly on compact subsets ofC\[—R, R], where Zi{x, s) is the solution to the algebraic equation 
h2.16\) with smallest modulus. 

Proof. The claim about the boundedness of the zeros of Qk,n follows in a rather standard way 
from the assumptions, see e.g. Proof of Lemma 2.2]. Now we turn to the claim (|4.ip . We 
consider the family of functions 

= (-^^^^ I fc,n€N,fc< (,s + l)nl. (4.2) 

yQk+l,n{^) J 

From the assumption that the zeros of {Qk,n)k are real and interlacing, it follows that "H is a 
normal family (in the sense of Montel) on C \ [—R, R], see e.g. [26, Proof of Lemma 2.2]. 
Using induction on I, we will show the following. 

Claim: For any m g {0, . . . , r — 1} and for each I > 1, the following holds. If {ki)i, {ni)i are 
sequences of non- negative integers with ki,ni -> oo, rki/rii s as i — > oo, so that 

f[x) := hm 



Qrki-\-m-\-r,ni (^) 

exists for |a;| > R, then 

f{x)^z,{x,s){l + Oix-')) 

as a; — )■ oo. 

Let us prove this claim. We have zi{x, s) = x~^{l + 0{l/x)) as x —>■ oo fLemma l2.3f a)). and 
so it is clear that the claim holds for ^ = 1. 

Now assume that the claim holds for Z > 1. We will prove that it also holds for I + r. We 
will prove this when to = 0; the proof for the other values of m can be given in a similar way. 
Letting {ki)i, {ni)i be as in the claim, our goal will be to prove that the function e{x) defined by 

lim =.i(x,.)(l + 6(:r))-\ (4.3) 

* ^ rki-\-r ^Uiy^ } 

satisfies e(x) = ©(x"'^'') for x — oo. 

Let us prove this. Since rki/ui — > s as i cx), we may assume that 

rh < (s + l)n, -(r-1) 

for every i. For j = — (r — 1), . . . , r/3, we then have that 

Qrki —j,ni 
Qrki+l— j ^Ui 

belongs to the family %. Since "H is a normal family, we may assume, by passing to a subsequence 
if necessary, that 

/(■'•)(X)^ lim Qrk.-3,nAx) 
i^oo Llrki + l~j,ni\X) 
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exists for a: G C \ [—R, R] and j — — {r ~ 1), . . . , r(3. 

Taking the kith block row in (|2.10p with n — Ui and using (I2.1ip . we obtain the matrix- vector 
relation 



(4.4) 



rxr{l3+2) 



r(/3+2)xl 



where we denote with Qfe,n(a;) the r x 1 column vector 

(Qrk,n 
\ 
Qr{k+l)-l.n{x)j 

Dividing both sides of (j4.4p by the scalar function Qrki,ni{x), and taking the limit i — > oo, we 
find 

/Q(«(x) 



where we used (|2.13p and where we set 



Q(^)(a;) := lim 



Q/ci —j 



W-^Hx 



i^oc Qrki,ni{x) 



J 



-1,. ..,/?, 



(4.5) 



(4.6) 



which exists entrywise due to our assumptions. It will be convenient to rewrite Q^^-'(x), j 
0, . . . , /3, as a telescoping product: 



J-yoo Qrki,ni(X) 



,. Qk,-j,7ii{x) ,. Qfci-j + l.„, (x) Qfei,n, (x) 

hm — — / , Imi „ — , ' X • • • X hm ■ 



Qk,-j + l,nAx) Qfei-j+2,n. (a;) 



(4.7) 



where by abuse of notation we write ^ and ab for two vectors a, b of length r to denote their 
entrywise quotient and product respectively. Each of the limits in (j4.7|) exists again entrywise due 
to our assumptions. Applying the induction hypothesis to each of the limits in the telescoping 
product (|4.7p . we find that 



(x) = Zi{x, s (x), X ^ oo 



iO)( 



(4- 



for any j = 0, . . . , /?. 

Next we turn to the term Ai ^•'Q^~^^(x) in the expansion of the right hand side of (j4.5l) . 
Denoting ei := (1,0,..., 0)^, we can write this as 



^(-i)q(-i)(^) =^(-1) 



Qk, 



+ l,n. 



(x) 



— lim 



Qr(ki + \),ni (a;) 



yi^oo Qrk,,n,{x) 

= zi{x,s)-\l + e{x))Ai-~^'>ei 

= Zi(a:,s)-i(l + 6(x))4-i)QW(x), 



(4.9) 
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where the first step follows by definition, the second step follows since the matrix Ai is zero 
except for its bottom left entry, cf. (j2.12p . the third step is a consequence of (|4.3[) . and the last 
step uses that Q^°)(a;) in (j4.6p has its first entry equal to 1. 

Inserting (j4.8p - (14.9p in (|4.5p yields the matrix-vector relation 

2;Q(°)(a;) = S(a;,s)Q(°)(2^) (4-10) 

where the matrix B(x, s) satisfies 

B{x, s) = sf{l + 0{x-')) + ■■■ + Ai^h^ix, + 0{x-')) + 

+Ai-^hi{x,s)-\l + e{x))Y x-^oo. (4.11) 

We can rewrite (|4.1ip as 

B{x, s) = [Ail^'>zi{x, sf + --- + A^^hiix, s) + 4°) + Ai-^hi{x, 

+ Air^'>zi{x,s)-'^e{x) + 0{x-'-''), x-^oo, (4.12) 

by using that zi{x,s) = 0{x~'') as a; — > oo fLemma l2.3r a)). 
Relation (j4.10p clearly implies that 

det{Bix,s)-xIr) = 0. (4.13) 

Expanding the determinant (|4.13p for |a:| large, with the help of (j4.12p . we obtain 

det (^Ai^hiix, s f + ■■■ + Ai^hi{x, s) + + Ai-^hiix, s)-'^ - x/^) 

+ {-lY+h(x)x''{l + O{l/x)) + O{x-'):^0, x-^oo. (4.14) 

Here the terms in the second line of (14.141) can be justified by using the special structure of A^f^ 
and Ai in (|2.12p . and using again the fact that zi{x, s) = a;~''(l + 0{l/x)) as x — cx). 

Now the determinant in the first line of (|4.14p vanishes identically, since by definition zi = 
zi{x,s) is a root of (I^TTC)) . cf. So reduces to 

(-l)''+ie(a;)x'^(l + 0{l/x)) + ©(a;"') = 0, x ^ oo, 

which implies in turn that e{x) = ©(a;^'^'') as a; — > oo. This proves the induction step, thereby 
establishing the claim. 

Having proved the claim, the proof of Lemma |4. II can now be finished from a standard normal 
family argument as in [26J. □ 

With Lemma H?T] in place, the proof of Theorem 12.41 can be finished as in [2(3, Proof of Theo- 
rem 1.2]. □ 

Remark 4.2. The above proof also shows that 

^ Qrki ,ni (•^) 



Q(°)(x) lim ^ 



i^oo Qrki,niix) 



r{ki-\-l) — l.ni 
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satisfies 

a;Q(")(a;) ^ As{zi{x))Q^°\x), xeC\[-R,R], 

by virtue of ((27T4| and the fact that (|47T0| - (l4TT|) hold with I arbitrarily large. So Q(°)(a;) is a 
vector with first component equal to 1 which lies in the null space of the matrix As{zi{x)) — xir- 
In fact, it can be shown that there is a unique vector v(a;) satisfying this condition, for all 
X £C \ [—R, R], and hence we have 



QW(x)= lim -\- 

rk/n^s <^rk,n[X) 



^ Qrk.nix) 
\Qr{k+l)-l,n{x) J 



= v(a;), a; e C\ 



5 Proof of Theorem 12.6 



5.1 The sets ri(s) and T2{s) 

Recall the functions Zj ~ Zj{x, s), defined as the solutions to the algebraic equation (|2.25p ordered 
by increasing modulus, (|2.26p . and the definitions (|2.19p and (|2.27p of the sets ri(s) and r2(s), 
whose structure we now describe. 

Proposition 5.1. ri(s) = [~cis, cis] andT2{s) — iM.\(—ic2S,ic2s), where ci and C2 are explicit 
constants given by i2.28\) and \2.29\) . respectively. 

Proof. By (|2.25p . specializing (|2.16p to the present setting gives 

fs{z,x) = x^ - P{z,s), 

where 

, , U + zs^(l + P)fU + zs%^) , , 

P(z,s) = -^-^ ^ ^ ^, where6 = tant. (5.1) 

It is clear that x e ^j{s) if and only if —x e rj(s), so it will be convenient to consider for a 
moment the sets 

T'^{s):={y\y^x\xeT,is)}. (5.2) 

We begin by establishing that (s)U F^s) C M; the proof of this fact will be very similar 
to the proof of Lemma 4.1 in [26\. Suppose y S Fj(s) U ^2(3). We can assume without loss of 
generality that y is not a branch point, since the number of branch points is finite and Fi(s) and 
F2 (s) in our case have no isolated points [T71 . Thus there exist distinct zi,Z2 ^ C such that 
l^il = 1-^21 =: ''' and P{zi, s) = y = P{z2, s). By the factorization of P, we see that z 1— )■ P{z, s) 
has only negative real zeros and therefore the even function [— tt, tt] B 9 |P(re**, s)\ is strictly 
decreasing on (0,7r), which implies in turn that zi — '22. But 



y = Piz2,s) = P{zi,s) = P(zi, s) = y, 

so y is real, and hence F^(s)U F2(s) C M. This argument also shows that F^(s) n F|(s) may 
contain only branch points, since otherwise there would be three distinct values zi,Z2,Z3 with 
the same modulus and y = P{zi, s) = P{z2, s) = P{z3, s) which is clearly impossible. 

Studying the function M 9 z 1— )■ P{z, s) for fixed s > 0, we see that it has two local minima 
at the points 




22 



P{z,s) 

\ yiis)' 












ziis) ^ 

'y2{s) 





Figure 4: Plot of P{z, s), z e 



and 



with 




2 , 2754 + 1862 - 1 + \/6^TT(962 + 1)3/2 



3262 



(5.3) 



(cis)2 > 0, (5.4) 



and 



2/2(5) :=F(z2,s) = s' 



276^ + 1862 - 1 - ^/62^+T(962 + 1)^/2 
3262 



= (ic2s)2 < 0, 



(5.5) 



and that is a local maximum attained at the point 

4 



20 (s) 



1 + 62 



(see Fig. [3]). 
Since 



P{z,s) ^ y 



(5.6) 



is a polynomial equation in z with real coefficients, it has two complex conjugate solutions if 
y e (-00,2/2(3)) U (0,yi(s)). If y e R \ ((-00, ?/2(s)] U [0,yi(.s)]), then all three solutions to (lEB)) 
are real and it is easy to see that there can be at most a finite number of such y for which two 
of these roots have the same modulus. But as already mentioned, ri(s) and r2(s) cannot have 
isolated points, hence r2(s) ur2(s) C (—00,2/2(5)] U [0,2/i(s)]. 

Now consider the interval [0,2/i(s)]. For the branch point y — yi{s), we have that 2:1(5) > 
is a double root to (|5.6p . and there is also a negative root z^{s) which is smaller than 22(5) (see 
Fig. HI). Therefore 

\zi{s)\ - |z_(s)| < zi(s) + Z2{S) = < 0, 
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i.e., the real negative root z_(s) has larger modulus than the double positive root zi{s), so 
2/1 (s) G rf (s) \ r| (s) . In a similar way it follows that G (s) \ (s) and 2/2 (s) € r| (s) \ (s) . 
Thus the branch points are not in r^(s) n ^2(3), and therefore 

r?(s)nr2(s) = 0, 

since we already proved that non-branch points are not in the intersection. 
We have now established that 

Tl{s) U Tlis) = [0, (cis)2] u (-c^, -{C2sf] 

with 0, (cis)2 G rf(s) and ~ic2sf G r^(s). We also established that r^(s)nr^(s) = 0, and then 
it follows by continuity that 

r?(s) = [o,(ci.)2], rl{s) = {-oo,-ic2sn 

This completes the proof. 

□ 

5.2 Equilibrium problem 

Given a measure /i in the complex plane, define the logarithmic potential of ^, 

J \y-x\ 



For j = 1,2, define the measures 



11/ Zj+ix,s) z'j^{x,s) 



d^i'Jx) = -— ' - dx, (5.7) 

^ 2 27^^ \^Zj+(.T,s) z^_(a;,s)y ^ ^ 

where dx is the (complex) line element on Tj{s). 

Proposition 5.2. Fix s > 0. The pair (^f , /i|) is i/ie unique minimizer of the energy functional 

J in, v) = I{p., /i) + v) - I{p., v) 

among all pairs of positive measures such that supp/i C ri(s), /d/i = 1; supp C r2(s) 

and J di' — 1/2. 

For all X G C, (/if,/^!) satisfies the Euler-Lagrange variational conditions 



2W'^i (x)-U''Hx)-l' = log , 

\zi{x,s)\ 

-w^i (x) + m^^i (x) = -i log -l^'^^' "^^l 



(5. 



|Z2(X,S)| ' 

for some constant P. 

Proof See [17J. □ 
We will now integrate (j5.8p to get an equilibrium problem for (i^i, 2^2), where 

i/2 := / M2ds, (5.9) 
in analogy with the definition of vi in (|2.23p . 
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Proposition 5.3. For all complex x, the vector {^1,^2) of measures satisfies the following con- 
ditions. Firstly, 

2U''^{x)-W^{x)-l + V{x)>Q, (5.10) 

where 

and I is some constant, with equality in (jS.lOp if and only if x ^ [— Ci,ci]. Secondly, 

-W^x) +2U''^{x) >0, (5.11) 

with equality if and only j/x e iR\ {—ic2,ic2). 
Furthermore, < f where 

/•OO 

/i^ds. (5.12) 

Proof. By Proposition 15.11 the supports of the measures /if and are subsets of the real and 
imaginary hues respectively, which increase/decrease linearly in s. This means that we can 
integrate the logarithmic potentials of /if and /i| with respect to s and change the order of 
integration to obtain 

fl pi pcis 



f U''Hx)ds^ f r' log^^d/if(2/)ds 

Jo Jq J -CIS \y-x\ 



ci f\y\/ci 1 

/ log, -dsd^^l{y)^:U-'-{x) 

-ci Jo \y ~ ^\ 

and similarly for i'2- The integrated variational conditions (jS.Sp thus become 

2U''^{x)-U-'^{x)- ri^ds = -l /'log^^^ds, (5.13) 
Jo ^ Jo \zi{x,s)\ 

-W^{x) + '2U'^-{x) = -\ /'log^^^ds. (5.14) 
Furthermore, since by definition |z2(a^, s)| > |zi(a;,s)|, 

1 f\ \Z2{X,S)\ ^ ^ 1 \^2{X,S)\ ^ 

' log - — —ds < - / log - — r-j-ds — V(x) 



2 Jq \zi{x,s)\ 2 Jq \zi{x,s)\ 
with equality if and only if a; G ri(l) ~ [— ci, ci]. Clearly, 1^2 must satisfy 



\dx\ Jo \dt\'^'-Jo Wl'^'^W ^^-^^^ 

since x ^ supp/ij for s > \x\/c2. with equality in (j5.15p if and only if x G r2(l). Inserting into 
([gl^ and ((ET1) and putting I := I'ds gives the stated inequalities. 

□ 

Equations (|5.10p and (jS.lip are the Euler-Lagrange variational conditions for the equilibrium 
problem in Theorem 12.61 The external field V, density of the measure vi and upper constraint 
measure cr can be calculated explicitly, and the following subsections are devoted to these com- 
putations. 
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5.3 Calculation of V 

In the calculations that follow we make use of the function 
which is such that 

g(z) = lp(z/s2;s) (5.17) 
for every s > 0. So in particular Q{z) — P{z; 1). We also define 

Zj{x) := Zj{x, 1), j = 1,2,3, 
and these are the solutions of Q{z) — . Because of (j5.17p we have 

Zj{x, s) — —Zj , s > 0, (5.18) 

and it follows that 

= .>0, , = 1,2,3. (5.19) 

Zj[X,S) ox SZj[X/S) 

Proposition 5.4. The external field V is given by 

y(a;) = i /°°logM^i£}lds = (^-2<)|a:|, xER. (5.20) 
2 Jo \zi[x,s)\ 

Proof. Because of symmetry, we may assume a; > 0. 

It follows from Proposition 15.11 that Ti{s) is increasing with s and x € Ti{s) if and only if 
s > x/ci. Thus 1^2(3^, s)| = |zi(x, s)| for s > x/ci, and the integral that defines V{x) can be 
restricted to an integral over s G [0,x/ci]. Using (|5.19|) we obtain 



dV{x) 1 / 1 dz2{x,s) 1 dzi{x,s) 



dx 2 Jf^ \Z2{x,s) dx zi(x, s) dx 
1 r'^- ( ~z'2{x/s) z[{x/s) 



ds 



2 Jq \SZ2{x/ s) szi(x/s) 

1 r f -z'^iu) ~z[{u) 



ds 



2 \UZ2{u) uzi{u) 



du (5.21) 



where in the last step we made the change of variables u — x/s (recall that x > 0). Note that 
(|5.2ip does not depend on x. To evaluate (I5.2ip we note that both u 1— > zi{u) and u H> Z2{u) 
are one-to-one for u e [ci, 00) and they map the interval [ci, 00) onto (0, ii(ci)] and [z2(ci), 00), 
respectively. We split the integral (|5.2ip into two integrals, and apply a change of variables z = 
Zj(u) to each of them. Then combining the two integrals again, and noting that zi(ci) = Z2(ci) 
and that u — ^ Q{z) if z = Zj{u) with j = 1,2, we obtain 

dV(x) 1 f°° dz , , 

^ ' ' x>0. (5.22) 



(5.23) 



dx 2 Jo z^fW)' 
The integral in (I5.22p can be calculated explicitly, since 

4 + (62-i)z^ 



_d 

'dz 



arctan 



2^z(4 + 62z) 



(4+(l + 62)2)^z(4 + fe2^) 2z.jQ{z)' 
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Therefore 

dV(x) /4+(62-l)z\ f 4+(b^-l)z^ 

— lim arctan — — lim arctan — 

dx ^^0+ \2y'z{A + b^z) J ^^+°^ \_2^z(4 + 62z)^ 

= arctan 

2 \ 2b 

Using that b — tan t and applying trigonometric identities, we finally obtain 

dV^ 

— = 7r-2i. (5.24) 
da; 

Noting that V{0) = since G ^i{s) for every s > 0, we find the claimed expression for the 
external field by integrating (|5.24p with respect to x. □ 

Remark 5.5. Note that the external field V has the form to be expected by analyzing directly 
the asyniptotics of the weight functions Wi'. For the rescaled polynomials Qk,n the orthogonality 
conditions read 

/>oo 

Qk,n{x)x"^Wj{nx) dx — 0, m — 0, . . . ,kj — I, j — 1,2. 
Thus we have new effective weights Wjix) — Wj{nx). Using Stirling's formula, 



1 

2^^ 



_ ^n{2tjX-7r\x\){l+o{l)) 

Asymptotically, the dominant weight determining the potential associated with the distribution 
of zeros, will be max{wi{x),W2{x)} = e^"'-'^"^*^!^"^"'^^"^^''-', giving the external field V{x) — 
(tt - 2t) \x\. 

5.4 Density of ui 

Next we turn to the density of vi . 

Proposition 5.6. The measure vi is absolutely continuous with density given by (|2.30p - (j2.32p . 
Proof. By and (OTT)) . 

di^i 1 1 (z'-^^{x,s) z[_{x,sy 

^ ds, X € [— ci, ci . 



da; 2 27ri Ji^i/g^ \2;i+(a;, s) 2;i_(a;, s) 

We now make essentially the same calculations as in the proof of Proposition 15.41 Assuming 
X £ (0, ci) and using (j5.19p we obtain as in (|5.2ip 

da; 2 2TTtJ^ \uzi+(u) uzi^[u) J 

A change of variables z — zi±{u) leads to 

da; 2 27ri Ji^^(^^) z^Q{z) 
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which is an integral in the complex z-plane. Since x € (0, ci), we have that zi^{x) and zi^{x) 
are each others complex conjugate, and it can be shown that Im5i+(a;) < 0. The integral in 
(|5.25p is from zi+(a;) in the lower half plane to its complex conjugate in the upper half plane 
along a path in C \ (— oo,0]. The square root ^ Q{z) is defined and analytic in C \ (— 46~^,0] 
and it is positive for z real and positive. 



Let 



4+ (6^ - l)z 
2.yz{4 + b^zT' 



z = zi^{x), < X < Cl. 



(5.26) 



Then by ^^^^ and 



dx 



(x) — f arctan(w(x)) — arctan(u)(a;)) 



The arctangent is understood here as an analytic function 



arctan w — — log — 
2 \ 1 



w e C \ ((—ioo, —i) U (i, ioo)). 



Thus 



dui 
dx 



ix) 



Atti ^ 

ZTT 



/ 1 — iw{x) 1 + iw{x) 
\l + iw{x) l-iw{x)^ 
1 + iw{x) 



1 — iw{x) 



Using the defining relation Q{zi^{x)) — x^ in (|5.26p we can see that (|5.26l) is equal to (|2.31[) - 
(|2.32p and then (I2.30p follows. Note also that ^w{x) < since vi has a positive density. □ 



5.5 Calculation of upper constraint measure 



Proposition 5.7. The upper constraint measure a in (j5.12p is a multiple of the Lebesgue measure 
on the imaginary axis with density 2t/TT. 

Proof. Because of symmetry it is enough to consider the density on the positive imaginary axis. 
Let X G iR+. Then 



da da 1 /-l-l/^^ fz'^+{x,s) z'^^{x,s) 



\dx\ 



dx 



Z2+{x,s) Z2~ix,s) 



ds 



and by a calculation as in the proof of Proposition 15.41 this is equal to 



\dx\ An 



\UZ2+{u) UZ2-{u) 



du. 



Now we are going to make the change of variables z = Z2±{u). Since Z2+(ioo) — —ioo, 
Z2-(ioo) = +ioo, and Z2+{ic2) ~ Z2-{ic2), we obtain 



da 



(x) 



1 

An 



-\-ioo 



dz 



with integration along a contour from —ioo to 



o z{Q{z)y/^ 

ioo that intersects the real line in 



(5.27) 



1 



Z2±{iC2) = 22(1) = I 1 



962 + I 

52 + 1 
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and Figure|3]for ^2(5). It is easy to check (and this can also be seen from Figure^]) that 



62 I ' \l P + l J 62 + 1 ■ 

and so the contour intersects the real line in a point lying in between the simple root and the 
double root of Q{z). The branch of ((3(z))^/2 that is used in (|5.27p is the one that is defined and 
analytic in C \ (—00, —4/62] y jg^ ^^-j ^^^^^ ^^^^^ jg ^^j. ^ _ Z2±{ic2). Thus 

(Q(z))V2 ^ -^(4 + (1 + 62).)li±^^ (5.28) 

with principal branches of the fractional powers. 

We now deform the contour in (|5.27p to the positive real axis. In doing so, we pick up a 
residue contribution at z = — pqrj, which is 

-^(27rz)Res( — -i-— ,z = -— ^) =1, (5.29) 

since the residue turns out to be —2i, and so is independent of 6. The deformed contour goes 
along the positive real axis, starting at +00 on the lower side, and ending at +00 on the upper 
side of the real axis. It gives the contribution 

1 /■+°° / 1 1 \ . 



With our choice (|5.28p of square root we have {Q{z))]l'^ = —{Q{z))^_l'^ = ^ Q{z) where is the 
positive square root of a real and positive number. So the contribution from the positive real 
line is the integral 

1 /■+°° dz 



27T Jo Z^fW) 

and this integral we already calculated in (|5.22|) and (j5.24p . Its value is 

1 /■+°° dz 2t 



1 + — . (5.30) 



2^ 7o z^fW) 
Adding ((Og)) and ((OO)) we find 

da _ 2t 

\dx\ TT 

and this proves the proposition. □ 
5.6 Proof of Theorem 12.61 

The equilibrium problem has a unique minimizer which will satisfy the Euler- Lagrange variational 
conditions, and Proposition 15.31 shows that (i^i,i^2) has this property. The explicit expressions 
for the external field, upper constraint measure and density of vi are calculated in Propositions 
OOandEi □ 
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